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Dalitz plot analysis is an excellent way to study the 
dynamics of three-body charm decays. These decays are 
expected to proceed predominantly through intermediate 
quasi-two-body modes |l[ and experimentally this is the 
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observed pattern. Dalitz plot analyses can provide new 
information on the resonances that contribute to the ob- 
served three-body final states. In addition, since the in- 
termediate quasi-two-body modes are dominated by light 
quark meson resonances, new information on light meson 
spectroscopy can be obtained. 

Some puzzles still remain in light meson spectroscopy. 
There are new claims for the existence of broad states 
close to threshold such as k(800) and / (600) [2]. The 
new evidence has reopened discussion of the composition 
of the ground state J PC = ++ nonet, and of the possi- 
bility that states such as the a (980) or / (980) may be 
4-quark states, due to their proximity to the KK thresh- 
old Q. This hypothesis can be tested only through ac- 
curate measurements of the branching fractions and the 
couplings to different final states. It is therefore impor- 
tant to have precise information on the structure of the 
7T7T 5-wave. In addition, comparison between the pro- 
duction of these states in decays of differently flavored 
charmed mesons D°(cu), D + (cd) and Df(cs) can yield 
new information on their possible quark composition. In 
this context, Df mesons can give information on the 
structure of the scalar amplitude coupled to ss. Another 
benefit of studying charm decays is that, in some cases, 
partial wave analyses are able to isolate the scalar con- 
tribution with almost no background. 

This paper focuses on the study of the three-body 
Df meson decays to ir + ir~ir + [4| and performs, for the 
first time, a Model-Independent Partial Wave Analysis 
(MIPWA) [5(. Previous Dalitz plot analyses of this de- 
cay mode were based on much smaller data samples d, 0] 
and did not have sufficient statistics to perform the de- 
tailed analysis reported here. 

This paper is organized as follows. Section II briefly 
describes the BABAR detector, while Section III gives de- 
tails on event reconstruction. Section IV is devoted to 
the evaluation of the selection efficiency. Section V deals 
with the Dalitz plot analysis of Df — > ir + ir~ir + and re- 
sults are given in Section VI. The measurement of the 
Df branching fraction is described in Section VII. 



II. THE BABAR DETECTOR AND DATASET 

Tha analysis is based on data collected with the BABAR 
detector at the PEP-II asymmetric-energy e + e~ collider 
at SLAC . The data sample used in this analysis corre- 
sponds to an integrated luminosity of 347.5 fb" 1 recorded 
at the T(45) resonance (on-peak) and 36.5 fb _1 collected 
40 MeV below the resonance (off-peak) . The BABAR de- 
tector is described in detail elsewhere Q. The follow- 
ing is a brief summary of the components important to 
this analysis. Charged particles are detected and their 
momenta measured by a combination of a cylindrical 
drift chamber (DCH) and a silicon vertex tracker (SVT), 
both operating within a 1.5 T solenoidal magnetic field. 
A ring-imaging Cherenkov detector (DIRC) is used for 
charged-particle identification. Photon energies are mea- 



sured with a Csl electromagnetic calorimeter (EMC). In- 
formation from the DIRC and energy-loss measurements 
in the DCH and SVT are used to identify charged kaon 
and pion candidates. Monte Carlo (MC) events used in 
this analysis, e + e _ — > cc, are generated using the JET- 
SET program @, and the generated particles are prop- 
agated through a model of the BABAR detector with the 
GEANT4 simulation package [lOj. Radiative corrections 
for signal and background processes are simulated using 
PHOTOS [Hi. 

III. EVENT SELECTION AND Df 
RECONSTRUCTION 

Events corresponding to the three-body decay: 

Df -> TT+TT-TT+ (1) 

are reconstructed from the sample of events having at 
least three reconstructed charged tracks with a net charge 
of ±1 and having a minimum transverse momentum of 
0.1 GeV/c. Tracks from Df decays are identified as pions 
or kaons by the Cherenkov angle 8 C measured with the 
DIRC. The typical separation between pions and kaons 
varies from 8a at 2 GeV/c to 2.5<r at 4 GeV/c, where a is 
the average resolution on 9 C . Lower momentum kaons are 
identified with a combination of 9 C (for momenta down to 
0.7 GeV/c) and measurements of ionization energy loss 
dE/dx in the DCH and SVT. The particle identification 
efficiency is w 95%, while the misidentification rate for 
kaons is ks 5%. Photons are identified as EMC clusters 
that do not have a spatial match with a charged track, 
and that have a minimum energy of 100 MeV. To reject 
background, the lateral energy is required to be less than 
0.8. The three tracks are fitted to a common vertex, and 
the x 2 fit probability (labeled Pi) must be greater than 
0.1 %. A separate kinematic fit which makes use of the 
Df mass constraint, to be used in the Dalitz plot anal- 
ysis, is also performed. To help discriminate signal from 
background, an additional fit which uses the constraint 
that the three tracks originate from the e + e~ luminous 
region (beam spot) is performed. We label the \ 2 prob- 
ability of this fit as P2, and it is expected to be large 
for background and small for Df signal events, since in 
general the latter will have a measurable flight distance. 

The combinatorial background is reduced by requiring 
the Df to originate from the decay 

r>:(2112)+-£>+ 7 (2) 

using the mass difference Am = m(Tr + ir~n + "f) — 
m(7r + 7r~7r + ). We cannot reliably extract the Z?*(2112) + 
characteristics using the 3tt decay mode due to the large 
background below the signal peak. Therefore we use the 
decay 

Df -» K + K-ir + , (3) 

which has a much larger signal to background ratio. 
Fitting the mass difference Am = m(K + K~ir + j) — 
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m(K + K~ix + ) for this decay mode with a polynomial de- 
scribing the background and a single Gaussian for the 
signal, we obtain a width a — 5.51 ± 0.04 MeV/c 2 . Since 
the experimental resolution in Am is similar for the two 
-D+ decay modes, we require the value of Am for the 
D+ — > 7r + 7r~7r + mode to be within ±2cr of the Review 
of Particle Physics [12] value of the (D*(2112)+ - D+) 
mass difference. At this stage the three-pion invariant 
mass signal region, defined between (— 2a, 2a), where a 
is estimated by a Gaussian fit to the Df lineshapes, has 
a purity (signal/ (signal+background)) of 4.3%. 

Each D+ candidate is characterized by three variables: 
the center of mass momentum p* , the difference in proba- 
bility Pi — P2, and the signed decay distance d xy between 
the Df decay vertex and the beam spot projected in the 
plane normal to the beam collision axis. The distribu- 
tions for these variables for background are inferred from 
the — > 7T + 7r~7r + invariant mass sidebands defined 
between (— 9a, — 5a) and (5a, 9a). Since these variables 
are (to a good approximation) independent of the de- 
cay mode, the distributions for the three-pion invariant 
mass signal, are inferred from the Df — * K + K~ir + de- 
cay. These normalized distributions arc then combined 
in a likelihood ratio test. The cut on the likelihood ratio 
has been chosen in order to obtain the largest statistics 
with background small enough to perform a Dalitz plot 
analysis. 

Many possible background sources are examined. A 
small background contribution due to the decay D* + — > 
ir + D° where D° — > ir + ir~ is addressed by removing 
events with |to(7t + 7t - ) — m D o\ < 20.7 MeV/c 2 and 
TO(7r+7r-7r+) - to(7t+7t-) < 0.1475 GeV/c 2 . Particle 
misidentification, in which a kaon (K m i S ) is wrongly iden- 
tified as a pion, is tested by assigning the kaon mass to 
each pion in turn. In this way we observe a clean signal 
in the mass difference m(Tr + K~ is TT + ) — m(ir + K~ is ) due 
to the decay D*+ -» tt+D where D° -> K mis ir+ . Re- 
moving events with \m(K mis ir + ) — rn D o \ < 21.7 MeV/c 2 
and m(7r+if mis 7r+) - m(K mis TT+) < 0.1475 GeV/c 2 di- 
minishes this background. Finally, events having more 
than one candidate are removed from the sample (1.2 % 
of the events). 

The resulting 7r + 7r~7r + mass distribution is shown in 
Fig. [lja). This distribution has been fitted with a single 
Gaussian for the signal and a linear background function. 
The fit gives a D+ mass of (1968.1 ± 0.1) MeV/c 2 and 
width a = 7.77 ± 0.09 MeV/c 2 (statistical error only). 
The signal region contains 13179 events with a purity 
of 80%. The resulting Dalitz plot, symmetrized along 
the two axes, is shown in Fig. [Tf b) . For this distribu- 
tion, and in the following Dalitz plot analysis, we use the 
track momenta obtained from the Df mass-constrained 
fit. We observe a clear / (980) signal, evidenced by the 
two narrow crossing bands. We also observe a broad ac- 
cumulation of events in the 1.9 GeV 2 /c 4 region. 




2.05 



0.5 1 1.5 2 2.5 3 3.5 

m 2 (7i + 71") (GeV 2 /c 4 ) 

FIG. 1: (a) ir + n~n + invariant mass distribution for the Df 
analysis sample. The line is the result of the fit described in 
the text, (b) Symmetrized Df — > n + n~n + Dalitz plot (two 
entries per event). 



IV. EFFICIENCY 

The efficiency for this decay mode is determined 
from a sample of Monte Carlo events in which the Df de- 
cay is generated according to phase space (i.e. such that 
the Dalitz plot is uniformly populated) . These events are 
passed through a full detector simulation and subjected 
to the same reconstruction and event selection procedure 
applied to the data. The distribution of the selected 
events in the Dalitz plot is then used to determine the to- 
tal reconstruction and selection efficiency. The MC sam- 
ple e + e~ — > D* + X, where D'* + — > "fDf, used to compute 
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this efficiency consists of 27.4 x 10 6 generated events for 
D+ -» 7T+7r-7r+ and 4.2 xlO 6 for D+ -» K+R--K+. The 
Dalitz plot is divided into small cells and the efficiency 
distribution is fitted with a second-order polynomial in 
two dimensions. The efficiency is found to be almost uni- 
form as a function of the 7r + 7r~ invariant mass with an 
average value of rj 1.6 %. This low efficiency is mainly 
due to the likelihood ratio selection: it is 18.0% without 
this cut. 

The experimental resolution as a function of the 7r + 7r~ 
mass has been computed as the difference between MC 
generated and reconstructed mass. It increases from 1.0 
to 2.5 MeV/c 2 from the tt+tt" threshold to 1.0 GeV/c 2 . 



V. DALITZ PLOT ANALYSIS 

An unbinned maximum likelihood fit is performed on 
the distribution of events in the Dalitz plot to determine 
the relative amplitudes and phases of intermediate res- 
onant and nonresonant states. The likelihood function 
is: 



events 



+ (1 — x{m)) 



E t \h\ 2 B? 



(4) 



where: 

• m\ and m 2 , are the squared 7r + 7r - effective masses; 

• x(m) is the mass-dependent fraction of signal, de- 
fined as x(m) — Q(m)+p{ m ) • Here G(m) and P(m) 
represent the Gaussian and the linear function used 
to fit the 7r + 7r~7r + mass spectrum. 

• r\(m\,rn^) is the efficiency, parametrized with a 
two-dimensional second order polynomial. 

• Ai,Bi describe signal and background amplitude 
contributions respectively. 

• ki are real factors describing the structure of back- 
ground. They are computed by fitting the sideband 
regions. 

• iAiA* — f AiA*ri{m\,m2)dm\drn\ and I B i are the 
normalization integrals for signal and background 
respectively. The products of efficiency and ampli- 
tudes are normalized using a numerical integration 
over the Dalitz plot. 

• Ci are complex coefficients allowed to vary during 
the fit procedure. 

The efficiency-corrected fraction due to the resonant 
or nonresonant contribution i is defined as follows: 



fi = 



\d\ 2 J \ Ai\ 2 dm\dm\ 
Ei,k c j c tl AjAldmidml 



The fi values do not necessarily add to 1 because of inter- 
ference effects. The uncertainty on each fi is evaluated 
by propagating the full covariance matrix obtained from 
the fit. 

The phase of each amplitude (i.e. the phase of 
the corresponding a) is measured with respect to the 
/2(1270)7r + amplitude. Each P-wave and 2?- wave am- 
plitude Ai is represented by the product of a complex 
Breit-Wigner (BW(m)) and a real angular term: 



A = BW{m) x T(fi). 



(6) 



where m is the 7r + 7r~ mass. The Breit-Wigner function 
includes the Blatt-Weisskopf form factors [13[ . The an- 
gular terms T(Q) are described in Ref. [Hj]. 

For the ir + ir~ 5-wave amplitude, we use a different 
approach because: 

• Scalar resonances have large uncertainties. In ad- 
dition, the existence of some states needs confirma- 
tion. 

• Modelling the 5-wave as a superposition of Breit- 
Wigners is unphysical since it leads to a violation 
of unitarity when broad resonances overlap. 

To overcome these problems, we use a Model- 
Independent Partial Wave Analysis introduced by the 
Fermilab E791 Collaboration [f|: instead of including 
the 5-wave amplitude as a superposition of relativistic 
Breit-Wigner functions, we divide the it + tt~ mass spec- 
trum into 29 slices and we parametrize the 5-wave by an 
interpolation between the 30 endpoints in the complex 
plane: 



A 



S — wave 



(m™) = Interp(c fe (r7i W7r )e^ fc(m '" r) ) fe= i 



,30- 



(7) 



(5) 



The amplitude and phase of each endpoint are free pa- 
rameters. The width of each slice is tuned to get ap- 
proximately the same number of ir + n~ combinations 
(~ 13179 x 2/29). Interpolation is implemented by a 
Relaxed Cubic Spline The phase is not constrained 
in a specific range in order to allow the spline to be a 
continuous function. 

The background shape is obtained by fitting the D+ 
sidebands. In this fit, resonances are assumed to be in- 
coherent, i.e. are represented by Breit-Wigner intensity 
terms only. A good representation of the background in- 
cludes contributions from Kg, p°(770) and three ad- hoc 
scalar resonances with free parameters. 

Resonances are included in sequence, keeping only 
those having a fractional significance greater than two 
standard deviations. 



VI. RESULTS 

The fit results (fractions and phases) are summarized 
in Table [T] The resulting 5-wave ir + ir~ amplitude and 
phase is shown in Fig.r2£a),(b) and is given numerically in 
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FIG. 2: (a) <S-wave amplitude extracted from the best fit, (b) corresponding iS-wave phase, (c) 5-wave amplitude compared to 
the FOCUS and E791 amplitudes, (d) 5-wave phase compared to the FOCUS and E791 phases. Errors are statistical only. 



Table HI1 Fig.[2fc),(d) show a comparison with the result- 
ing S- wave from the E791 experiment, which performed 
a Dalitz plot analysis using an isobar model [6[ , and the 
FOCUS experiment, which made use of the K-matrix for- 
malism [7j. In the two figures, the two bands have been 
obtained by propagating the measurement errors and as- 
suming no correlations. This assumption may influence 
the calculation of the uncertainties on the phases and 
amplitudes which are different in the two experiments. 



TABLE I: Results from the Df — > -k + tt~-k + Dalitz plot analy- 
sis. The table reports the fit fractions, amplitudes and phases. 
Errors are statistical and systematic respectively. 



Decay Mode Decay fraction(%) 


Amplitude 


Phase (rad) 


/ 2 (1270)7r+ 
p(770)tt+ 
p(1450)tt + 
<5>-wave 


10.1±1.5±1.1 
1.8±0.5±1.0 
2.3±0.8±1.7 

83.0±0.9±1.9 


1. (Fixed) 
0.19±0.02±0.12 
1.2±0.3±1.0 
Table |TT] 


0. (Fixed) 
1.1±0.1±0.2 
4.1±0.2±0.5 

Table [XT] 


Total 
X 2 /NDF 


97.2±3.7±3.8 

437 _ 1 r, 

422-64 







The Dalitz plot projections together with the fit results 
are shown in Fig. [3l Here we label with m 2 (n + ir~)\ ovl 
and ?n 2 (7r + Tr _ )high the lower and higher values of the 
two 7r + 7r~ mass combinations. 

The fit projections are obtained by generating a large 
number of phase space MC events [15| , weighting by the 
fit likelihood function, and normalizing the weighted sum 
to the observed number of events. There is good agree- 
ment between data and fit projections. Further tests of 
the fit quality are performed using unnormalized Y® mo- 
ment projections onto the 7r + 7r~ axis as functions of the 
helicity angle 9, which is defined as the angle between the 
7r~ and the in the tt + it~ rest frame (or tt + for D~) 
(two combinations per event). The ir + ir~ mass distribu- 
tion is then weighted by the spherical harmonic Y^(cos 6) 
(L = 1 — 6). The resulting distributions of the (Y°) 
are shown in Fig. 0J A straightforward interpretation 
of these distributions is difficult, due to reflections orig- 
inating from the symmetrization. However, the squares 
of the spin amplitudes appear in even moments, while 
interference terms appear in odd moments. 

The fit produces a good representation of the data 
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FIG. 3: Dalitz plot projections (points with error bars) and fit results (solid histogram), (a) m 2 (7r + 7r )i ow , (b) m 2 {-K + -K )hi g h, 
(c) total m 2 (7r + 7r~), (d) m 2 (n + 'K + ). The hatched histograms show the background distribution. 



for all projections. The fit \ 2 is computed by divid- 
ing the Dalitz plot into 30x30 cells with 422 cells hav- 
ing entries. We obtain x 2 / NDF = 437/(422 - 64) = 
1.2. The x 2 is also calculated using an adaptive bin- 
ning with an average number of events per cell ~35 
(X 2 /NDF = 365/(391 - 64) = 1.1), obtaining a X 2 prob- 
ability of 7.2%. 

Attempts to include other resonant contributions, such 
as w(782) or /2(1525), do not improve the fit quality. MC 
simulations have been performed in order to validate the 
method and test for possible multiple solutions. 

The results from the Dalitz plot analysis can be sum- 
marized as follows: 

• The decay is dominated by the £)+ — » 
(7r + 7r~)s_ wavo 7r + contribution. 

• The 5-wave shows, in both amplitude and phase, 
the expected behavior for the Jo (980) resonance. 

• The S-wave shows further activity, in both ampli- 
tude and phase, in the regions of the /o(1370) and 
/o(1500) resonances. 

• The 5-wave is small in the /o(600) region, indicat- 
ing that this resonance has a small coupling to ss. 



• There is an important contribution from Dj — > 
/2(1270)7r + whose size is in agreement with that 
reported by FOCUS, but a factor two smaller than 
that reported by E791. This is the largest contri- 
bution in charm decays from a spin-2 resonance. 

• We observe a similar trend for the 5-wave am- 
plitude and phase between the three experiments. 
Our results agree better (within uncertainties) with 
the results from FOCUS than those from E791. 

Our results may be compared with different measure- 
ments of the 7T7T amplitude and phase from many other 
sources. For a recent review, see (16j |. 

Systematic uncertainties on the fitted fractions are 
evaluated in different ways: 

• The background parametrization is performed us- 
ing the information from the lower/higher sideband 
only or both sidebands. 

• The Blatt-Weisskopf barrier factors have a single 
parameter r which we take to be 1.5 (GeV/c) -1 and 
which has been varied between and 3 (GeV/c) -1 . 

• Results from fits which give equivalent Dalitz plot 
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FIG. 4: Unnormalized spherical harmonic moments (Y°) as a function of tt + it effective mass. The data are presented with 
error bars, the histograms represent the fit projections. 



descriptions and similar sums of fractions (but 
worse likelihood) are considered. 

• The likelihood cut is relaxed but the mass cut on 
the 7r + 7r - 7r + is narrowed in order to obtain a sim- 
ilar purity. 

• The purity of the signal, the resonance parameters 
and the efficiency coefficients arc varied within their 
statistical errors. 

• The p(770) and p(1450) parametrization is modi- 
fied according to the Gounaris-Sakurai model [l7T |. 

• The number of steps used to describe the S-wave 
has been varied by ±2. 



VII. BRANCHING FRACTION 

Since the two Df decay channels (1) and (3) have simi- 
lar topologies, the ratio of branching fractions is expected 
to have a reduced systematic uncertainty. We therefore 
select events from the two decay modes using similar 
selection criteria for the Z?* + selection and for the likeli- 
hood test. For this measurement, a looser likelihood cut 
is used. 



The ratio of branching fractions is evaluated as: 



BR 



^x,y e (x,y) 



(8) 



where Ni{x, y) represents the number of events measured 
for channel i, and £i(x,y) is the corresponding efficiency 
in a given Dalitz plot cell (x, y). For this calculation each 
Dalitz plot was divided into 50x50 cells. 

To obtain the yields and measure the relative branch- 
ing fractions, the 7r + 7r~7r + and K + K~n + mass distri- 
butions are fit assuming a double Gaussian signal and 
linear background where all the parameters are floated. 
Systematic uncertainties, summarized in Table IIII1 take 
into account uncertainties from MC statistics and from 
the selection criteria used. 

The resulting ratio is: 



) 



0.199 ±0.004 ±0.009 



(9) 



B(D+ ^ K+K-TT+) 

consistent, within one standard deviation, with the 
PDG [12] value: 0.265 ± 0.041 ± 0.031. It is also con- 
sistent with a recent measurement from CLEO [l8j : 
0.202 ±0.011 ±0.009. 

The study of the Df — > K + K~n + decay can give 
new information on the KK 5-wave. This information 
together with the results reported in this analysis will 
enable new measurements of the /o(980) couplings to 
tttt/KK. 
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VIII. CONCLUSIONS 

A Dalitz plot analysis of approximately 13,000 Df — > 
7r + 7r _ 7r + has been performed. The fit measures frac- 
tions and phases for quasi- two-body decay modes. 
The amplitude and phase of the 7r + 7r~ 5-wave is ex- 
tracted in a model-independent way for the first time. 
We also measure with high precision the B(Df — > 
TT+TT-n+)/B(Df -► K+K-1T+) ratio. 

TABLE II: Amplitude and phase of the ir~ir + <S-wave ampli- 
tude determined with the MIPWA fit described in the text. 
The first error is statistical while the second is systematic. 



Interval Mass ( GeV/c 2 ) Amplitude Phase (radians) 



1 


0.28 


2.7 


± 


1.5 


± 


2.4 


-3.4 


± 


1.0 


± 


1.3 


2 


0.448 


2.2 


± 


1.2 


± 


1.3 


-3.9 


± 


0.5 


± 


0.4 


3 


0.55 


3.2 


± 


0.8 


± 


1.1 


-3.7 


± 


0.3 


± 


0.3 


4 


0.647 


3.3 


± 


0.7 


± 


0.9 


-3.7 


± 


0.2 


± 


0.3 


5 


0.736 


5.0 


± 


0.7 


± 


1.1 


-3.4 


± 


0.1 


± 


0.2 


6 


0.803 


5.1 


± 


0.7 


± 


0.8 


-2.9 


± 


0.1 


± 


0.2 


7 


0.873 


6.7 


± 


0.7 


± 


0.7 


-2.6 


± 


0.1 


± 


0.3 


8 


0.921 


10.7 


± 


1.0 


± 


0.9 


-2.2 


± 


0.1 


± 


0.2 


9 


0.951 


16.3 


± 


1.6 


± 


1.2 


-1.9 


± 


0.1 


± 


0.2 


10 


0.968 


22.9 


± 


2.3 


± 


1.5 


-1.4 


± 


0.1 


± 


0.1 


11 


0.981 


27.2 


± 


2.7 


± 


1.6 


-0.8 


± 


0.1 


± 


0.2 


12 


0.993 


20.4 


± 


2.0 


± 


0.9 


-0.3 


± 


0.1 


± 


0.2 


13 


1.024 


11.8 


± 


1.2 


± 


0.5 


0.1 


± 


0.1 


± 


0.2 


14 


1.078 


8.8 


± 


0.9 


± 


0.3 


0.4 


± 


0.1 


± 


0.1 


15 


1.135 


7.4 


± 


0.7 


± 


0.3 


0.9 


± 


0.1 


± 


0.1 


16 


1.193 


6.3 


± 


0.5 


± 


0.2 


1.1 


± 


0.1 


± 


0.1 


17 


1.235 


7.0 


± 


0.5 


± 


0.3 


1.4 


± 


0.1 


± 


0.1 


18 


1.267 


6.9 


± 


0.5 


± 


0.3 


1.4 


± 


0.1 


± 


0.1 


19 


1.297 


6.1 


± 


0.6 


± 


0.6 


1.8 


± 


0.1 


± 


0.1 


20 


1.323 


6.7 


± 


0.6 


± 


0.5 


1.7 


± 


0.1 


± 


0.1 


21 


1.35 


7.0 


± 


0.8 


± 


0.6 


1.8 


± 


0.1 


± 


0.2 


22 


1.376 


7.5 


± 


0.8 


± 


0.7 


2.0 


± 


0.1 


± 


0.1 


23 


1.402 


9.2 


± 


1.0 


± 


0.9 


2.1 


± 


0.1 


± 


0.1 


24 


1.427 


9.1 


± 


1.0 


± 


0.9 


2.3 


± 


0.1 


± 


0.2 


25 


1.455 


9.1 


± 


1.0 


± 


1.6 


2.6 


± 


0.1 


± 


0.1 


26 


1.492 


7.0 


± 


0.9 


± 


1.1 


3.1 


± 


0.1 


± 


0.2 


27 


1.557 


2.3 


± 


0.5 


± 


0.7 


4.3 


± 


0.2 


± 


0.4 


28 


1.64 


2.8 


± 


1.1 


± 


1.3 


4.7 


± 


0.3 


± 


0.7 


29 


1.735 


3.1 


± 


1.1 


± 


2.3 


6.0 


± 


0.5 


± 


1.4 



TABLE III: Summary of systematic uncertainties on the 
B{Dt -> Tv + Ti-Tv + )/B(Df -f K+K--K+) ratio. 



Source 


Systematic uncertainties (%) 


MC statistics 


0.9 


Am cut 


1.5 


Likelihood cut 


2.6 


Particle identification 


3.0 


Total 


4.3 



IX. ACKNOWLEDGEMENTS 



We are grateful for the extraordinary contributions of 
our PEP-II colleagues in achieving the excellent luminos- 
ity and machine conditions that have made this work pos- 
sible. The success of this project also relies critically on 
the expertise and dedication of the computing organiza- 
tions that support BABAR. The collaborating institutions 
wish to thank SLAC for its support and the kind hospi- 
tality extended to them. This work is supported by the 
US Department of Energy and National Science Foun- 
dation, the Natural Sciences and Engineering Research 
Council (Canada), the Commissariat a PEncrgic Atom- 
ique and Institut National de Physique Nucleaire et de 
Physique des Particules (France) , the Bundcsministcrium 
fur Bildung und Forschung and Deutsche Forschungsge- 
meinschaft (Germany), the Istituto Nazionale di Fisica 
Nuclearc (Italy), the Foundation for Fundamental Re- 
search on Matter (The Netherlands) , the Research Coun- 
cil of Norway, the Ministry of Education and Science of 
the Russian Federation, Ministerio de Education y Cien- 
cia (Spain), and the Science and Technology Facilities 
Council (United Kingdom). Individuals have received 
support from the Marie-Curie IEF program (European 
Union) and the A. P. Sloan Foundation. 



[1] M. Bauer, B. Stech and M. Wirbel, Z. Phys. C 34, 103 
(1987). 

[2] E. M. Aitala et al. [E791 Collaboration], Phys. Rev. Lett. 
89, 121801 (2002); E. M. Aitala et al. [E791 Collabora- 
tion], Phys. Rev. Lett. 86, 770 (2001). 

[3] See for example F. E. Close and N. A. Tornqvist, J. Phys. 
G 28, R249 (2002). 

[4] All references in this paper to an explicit decay mode 
imply the use of the charge conjugate decay also. 

[5] E. M. Aitala et al. [E791 Collaboration], Phys. Rev. D 
73, 032004 (2006). 



[6] E. M. Aitala et al. [E791 Collaboration], Phys. Rev. Lett. 

86, 765 (2001). 
[7] J. M. Link et al. [FOCUS Collaboration], Phys. Lett. B 

585, 200 (2004). 
[8] B. Aubert et al. [BABAR Collaboration], Nucl. Instrum. 

Meth. A 479, 1 (2002). 
[9] T. Sjostrand, S. Mrenna and P. Skands, JHEP 05, 026 

(2006). 

[10] S. Agostinelli et al. [GEANT4 Collaboration], Nucl. In- 
strum. Meth. A 506, 250 (2003). 
[11] P. Golonka et al., Comput. Phys. Commun. 174, 818 



13 



(2006). 

[12] W. M. Yao et al, Review of Particle Physics, J. Phys. G 
33, 1 (2006) and 2007 partial update for 2008. 

[13] J. M. Blatt and V. F. Weisskopf, Theoretical Nuclear 
Physics, John Wiley & Sons, New York, 1952. 

[14] K. S. Kolbig and H. Lipps, Cubic Splines and Their In- 
tegrals, CERN Program Library, E211. 

[15] F. James, N-Body Monte-Carlo Event Generator, CERN 



Program Library, W515. 
[16] E. Klempt and A. Zaitsev, Physics Reports 454 (2007). 
[17] G. J. Gounaris and J. J. Sakurai, Phys. Rev. Lett. 21, 

244 (1968). 

[18] J. P. Alexander et al. [CLEO Collaboration], Phys. Rev. 
Lett. 100, 161804 (2008). 



